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DYNAMIC  RESPONSE  OF  AN  ELASTIC  PLATE 
CONTAINING  PERIODIC  MASSES 


L  INTRODUCTION 


Plate  theory  has  been  researched  extensively  for  many  years.  Early  plate  and  beam  theory' 
modeled  displacement  in  thin  plates  and  beams.  These  early  models  contain  primarily  flexural 
wave  dynamics  and  are  inaccurate  at  high  frequencies  and  wavenumbers.  Rotary  inertia  and  shear 
effects^  were  added  to  the  flexural  wave  model  to  obtain  more  accurate  results  at  increased 
frequencies.  Fully  elastic  models^  were  developed  to  incorporate  plate  dynamics  as  the  wave¬ 
lengths  of  energy  that  propagate  in  the  plate  began  to  approach  and  surpass  the  thickness  of  the 
plate.  Analysis  has  also  been  conducted  on  the  dispersion  curves  of  these  systems,  particularly  in 
the  area  of  free  wave  propagation.”*'^  To  a  lesser  extent,  the  mode  shapes  of  these  systems  have 
been  studied  and  documented.'"’  " 

The  complexity  of  plate  models  has  increased  over  the  years  by  the  addition  of  stiffeners 
(ribs)  or  masses.  The  modeling  technique  for  adding  a  stiffener  to  the  plate  is  similar  to  that  for 
adding  a  discrete  mass  to  the  plate.  The  problem  of  a  fluid-loaded  infinite  thin  plate  with  infinite 
sets  of  parallel  stiffeners  excited  by  a  point  load  has  been  analyzed  in  a  study  that  modeled  the 
stiffeners  as  line  forces.'^  This  problem  was  extended  to  include  a  moment  exerted  by  the 
stiffeners  and  forcing  functions  of  plane  wave,  line,  and  point  forces.'^’  The  fluid-loaded 
infinite  plate  problem  was  reformulated  for  a  finite  number  of  equally  spaced  stiffeners'^  and 
was  further  studied  for  randomly  spaced  stiffeners.'^  The  problem  of  a  fluid-loaded,  aperiodic- 
stiffened  infinite  plate  has  also  been  addressed,'^  as  has  the  analysis  of  a  finite-sized  plate 
containing  concentrated  masses.'*  In  these  studies,'^  '*  the  plate  model  has  been  either  a  thin 
plate  or  a  thin  plate  with  rotary  inertia  and  shear  effects.  Finally,  the  problem  was  modeled  using 
finite  elements  to  produce  numerical  solutions.'" 

This  report  presents  an  analytical  model  that  incorporates  an  infinite  number  of  periodically 
spaced  masses  into  the  equations  of  elasticity  that  model  motion  and  stress  in  a  two-dimensional 
solid.  The  formulation  of  the  problem  begins  with  elasticity  theory  that  models  the  motion  in  the 
solid  as  a  combination  of  dilatational  and  shear  waves.  From  this  theory,  an  expression  for  plate 
displacement  is  obtained.  The  displacements  are  then  inserted  into  stress  relationships  that  are 
set  equal  to  the  forces  acting  on  the  structure  by  the  masses.  The  problem  is  then  written  as  a 
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dynamic  system,  in  matrix  form,  where  the  left-hand  terms  represent  the  zero-order  modes  and 
are  equal  to  an  infinite  number  of  right-hand  terms  that  represent  the  masses  acting  on  the 
structure  plus  a  term  that  models  the  plane  wave  forcing  function.  Rewriting  this  zero-order  term 
by  increasing  and  decreasing  the  mode  index  results  in  an  expression  for  the  higher-order  modes. 
The  integer  shift  property  is  then  applied  to  the  right-hand  side  of  all  of  the  terms,  resulting  in  an 
infinite  set  of  equations  that  model  the  wave  propagation  coefficients  of  all  the  modes  of  the 
structure.  This  set  of  equations  is  truncated  to  a  finite  number  of  terms,  and  a  solution  to  the 
displacement  and  stress  field  is  calculated.  Two  different  cases  are  examined:  (1)  where  the 
masses  are  on  the  edge  of  the  plate  and  (2)  where  the  masses  reside  within  the  interior  of  the 
plate.  The  solution  is  compared  to  a  previously  solved  problem  at  low  frequency  where  the 
wavelength  of  the  plane  wave  forcing  function  is  large  compared  to  the  thickness  of  the  plate. 

A  numerical  example  of  a  high-frequency  problem  is  included  and  discussed. 
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2.  ELASTIC  PLATE  WITH  MASSES  ALIGNED 
ON  THE  LOWER  SURFACE 


The  first  problem  analyzed  is  that  of  an  elastic  plate  with  discrete  (point)  masses  at  the 
bottom  edge,  as  shown  in  figure  1.  The  masses  on  the  bottom  of  the  plate  are  equally  spaced  a 
distance  of  L  (meters)  in  the  x-direction  and  each  has  a  mass  per  unit  length  of  M  (kg/m).  The 
plate  has  a  thickness  of  h  (meters)  and  is  loaded  on  the  top  surface  with  a  forcing  function.  The 
model  is  based  on  the  following  assumptions:  (1)  the  forcing  function  acting  on  the  plate  is  a 
plane  wave  at  any  definite  wavenumber  and  frequency,  (2)  the  corresponding  response  of  the 
plate  is  at  definite  wavenumber  and  frequency,  (3)  motion  is  normal  and  tangential  to  the  plate  in 
one  direction  (two-dimensional  system),  (4)  the  plate  has  infinite  spatial  extent  in  the  x-direction, 
(5)  the  masses  have  translational  degrees  of  freedom  in  the  x-  and  z-directions,  and  (6)  the 
particle  motion  is  linear. 


Figure  1.  Elastic  Plate  with  Periodic  Edge  Masses 


The  motion  of  the  elastic  plate  is  governed  by  the  equation^® 
+  (/l  +  /i)VV.u  =  /?— , 

dP- 


(1) 


where  p  is  the  density  (kg/m^),  X  and  p  are  the  Lame  constants  (N/m^),  /  is  time  (seconds), 

•  denotes  a  vector  dot  product,  and  u  is  the  two-dimensional  Cartesian  coordinate  displacement 
vector  of  the  plate.  Equation  (1)  can  be  expanded  and  rewritten  using  four  boundary  conditions. 
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The  normal  stress  at  the  top  of  the  plate  (z  =  b)  is  equal  to  the  opposite  of  the  forcing  function 
and  is  expressed  as 


du^{x,b,t)  duAx,b,t) 

r^^(x, b, /)  =  (A  +  2/y) ^  ’  =  -/(x, /) , 


dx 


(2) 


where  u^{x,z,t)  is  the  displacement  in  the  z-direction  (meters),  u^{x,z,t)  is  the  displacement  in 
the  x-direction  (meters),  and  f{x,t)  is  the  forcing  function  exciting  the  plate  expressed  in  force 
per  unit  area  (N/m  ).  The  tangential  stress  at  the  top  of  the  plate  is  zero  and  is  written  as 


T^(x,b,t)  =  /j 


du  {x,b,t)  du  {x,b,t) 


dz 


dx 


=  0 


(3) 


The  normal  stress  at  the  bottom  of  the  plate  (z  =  d)  due  to  the  forces  induced  by  the  motion  of 
the  masses  is  equal  to  the  summation  of  the  mass  multiplied  by  the  acceleration  in  the  z-direction 
times  the  spatial  delta  Dirac  function.  This  expression  is 


du,(x,a,t)  duAx.a.t)  d^u  (x,a,t) 

T^^ix,a,t)  =  (A  +  2^)  M - i^^Six-nL).  (4) 


dz 


dx 


n=—oo 


dd 


The  tangential  stress  at  the  bottom  of  the  plate  is  equal  to  the  summation  of  the  mass  multiplied 
by  the  acceleration  in  the  x-direction  times  the  spatial  delta  Dirac  function.  This  equation  is 
written  as 


T.(x,a,t)  =  M 


du  (x,a,t)  du^{x,a,t) 
— 1 - +  — l - 


dz 


dx 


d^uAx,a,t) 

=  Z  M - 

«=-oo  dr 


(5) 


It  is  noted  that  each  mass  can  also  be  attached  to  ground  with  a  parallel  spring  and  damper  by 
replacing  the  mass  times  vertical  acceleration  term 


M 


2 

d  u^(x,a,t) 

^^2 


(6a) 


with 

4 


(6b) 


d  u  {x,a,t)  du  {x, a, t) 

M - - +  P — - +  Ku^(x,a,t) 


dt 


in  equation  (4),  where  P  is  the  viscous  damping  coefficient  per  unit  length  (Ns/m^),  and  K  is  the 
spring  constant  per  unit  length  (NW).  This  system,  shown  in  figure  2,  corresponds  to  a 
periodically  damped  and  stiffened  plate. 


Figure  2.  Elastic  Plate  with  Periodic  Edge  Masses,  Dampers,  and  Springs 


Modeling  the  displacement  as  a  dilatational  wave  and  a  shear  wave  and  inserting  this  term 
into  equation  ( 1 )  results  in  the  plate  displacements 

u^{x,z,t)  =  U^{k,z,co)t\p{ikx)c\p{i(ot) 

=  \A{k,  co)ik  sin(«z)  +  B{k,  co)ik  cos(az) 

-  C(k,co)j3cos{fiz)  +  D{k,a})j3s\n(^z)]exp(ikx)exp(io)t) ,  (7) 


and 
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u^{x,z,t)  =  ?72(^,z,ty)exp(/fcc)exp(/6>0 

=  [A{k,co)a  cos(<2rz)  -  B{k,(o)asm{az) 

+  C{k,  co)ik  sin(;^z)  +  D{k,  co)ik  cos(y9z)]exp(/A3r)  exp(/<yO,  (8) 


where  A{k,co),  B{k,oS),  C{k,Q^),  and  D{k,a))  are  complex  wave  propagation  coefficients  of  the  plate 
(determined  by  solving  the  system  of  equations  below);  /  =  ;  co  is  the  frequency  (rad/sec); 

a  is  the  modified  wavenumber  (rad/m)  associated  with  the  dilatational  wave  and  is  expressed  as 


where  ^^is  the  dilatational  wavenumber  and  is  equal  to  co/c^j,  with  being  the  dilatational 
wavespeed  (m/s);  /?  is  the  modified  wavenumber  (rad/m)  associated  with  the  shear  wave  and  is 
expressed  as 

p  =  ^kl-k^  ,  (10) 

where  k^  is  the  shear  wavenumber  (rad/m)  equal  to  «/c^,  with  Cg  being  the  shear  wavespeed 
(m/s);  and  k  is  the  spatial  wavenumber  in  the  x-direction  (rad/m).  The  relationships  between  the 
wavespeeds  cj  and  and  the  Lame  constants  are  determined  by 


- 


A  +  2/i 


and 


(11) 


Equations  (2)  through  (5)  are  now  transferred  into  the  wavenumber-frequency  domain 
using  equations  (7),  (8),  and 

/(x,  t)  =  F  exp(/A:x:)exp(;fy/) ,  (13) 
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n=co 


Z 


n=-oo 


u  ^{x,a,t) 


S{x-nL)  = 


n=oo 


^  f/  ^  (nL,  a,  co)  t\p{-iknL) 


L«=-oo 


exp(/<yO,  (14) 


and 


n=oo 


^  d  It  (x,  a,  t) 

^  M - - S{x  -  nL)  -  -  Map' 


n  =  -oo 


n=oo 


^  U  ^  {nL,  a,  co)  exp(-/^wZ,) 


./7  =  -00 


exp(/(y/).  (15) 


Additionally,  Poisson’s  summation  formula  is  applied  to  the  right-hand  side  of  equations  (14) 
and  ( 1 5),  yielding  the  stresses  in  the  wavenumber-frequency  domain  as 


{k,b,co)  dU Ak,b,(o) 

\,{k,b,co)  =  {A  +  2ju) - ^  +  A - ^  =-F. 


dz 


dx 


(16) 


\^{k,b,co)  =  14 


^U(k,b,co)  dU(k,b,co) 


dz 


dx 


=  0, 


(17) 


\^{k,a,co)  =  {A  +  2^) 


dU.{k,a,co)  dU  (k,a,co) 


dz 


+  A- 


dx 


-Mco^  jj  ,,  2.7tn  . 

2^  U^{k  + - ,a,(o). 


«  =  -00 


and 


Xj^{k,a,o})  =  n 


dU  {k, a, co)  dU  {k, a, co) 


dz 


—  Mco^  27in 

-  >  U^{k  + - ,a,co). 

I  I 

«=-o0  ^ 


dx 


(18) 


(19) 


Equations  (7)  and  (8)  could  be  inserted  into  equations  (16)  through  (19);  however,  the 
substitution  alone  would  not  produce  a  solution  to  the  problem,  because  the  result  would  be 
divergent  expressions.  Instead,  equations  (7)  and  (8)  are  substituted  into  equations  (16)  through 
( 1 9)  and  rewritten  as 
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(20) 


with 


Inn, 


•••  )]  [U^^^yt)]  [U(*^(yt  +  — )] 

L  L 


(21) 


where  [A^°^(^)]  is  a  4  x  4  matrix  that  models  the  dynamics  of  the  structure  for  «  =  0,  x^°^  is  the  4  x  1 
vector  of  wave  propagation  coefficients  for  «  =  0,  +  2;Tn/Z-)]  is  the  4  x  oo  block-partitioned 

matrix  that  represents  the  periodic  mass  loading  on  the  structure  for  n  =  -oo  to  /?  =  oo,  x*-®^  is  the 
00  X  1  vector  of  wave  propagation  coefficients  for  «  =  -oo  to  «  =  oo,  and  f  is  the  4  x  1  vector  that 
models  the  plane  wave  excitation.  (The  entries  of  the  matrixes  and  vectors  in  equations  (20)  and 
(21)  are  listed  in  the  appendix.)  To  facilitate  a  solution  to  the  problem,  the  integer  shift  property 
of  the  infinite  summation  in  equation  (20)  is  used.  Equation  (20)  is  w-indexed  and  becomes 


[  {k  +  — )]x('")  =  +  —  +  ^)]x(±«>)  +  f 

L  L  L 

=  [U(**’^(A: -I- ^^)]x(*‘*^  -I- f . 

L 


(22) 


Once  the  [A^""*]  matrix  is  integer-indexed  and  the  displacement  load  matrix  indices  have  been 
shifted,  the  system  equations  can  be  rewritten  using  all  the  ^-indexed  modes  as 

A  =  U  +  F  ,  (23) 

where  A  is  a  block-diagonal  matrix  and  is  equal  to 
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0 


A  = 


0  [A<“V)] 


0 


(24) 


U  is  a  rank  deficient,  block-partitioned  matrix  and  is  written  as 


[U(  ^){k-—)]  [U^^^A:)]  [u(*\A:  +  — )] 


U  = 


[U'  [U^^^A:)]  [U^‘^(A:  +  ^)] 

JLj  Li 

[U^"^\A:-— )]  [U^^^A:)]  [u(^\A:  +  — )] 


(25) 


and  F  is  the  plane  wave  load  vector 


(26) 


The  0  term  in  equation  (24)  is  a  4  x  4  matrix  whose  entries  are  all  zeros.  Once  equation 
(23)  is  assembled,  the  wave-propagation  coefficients  that  reside  in  the  x  vector  can  be 
determined  using 

=[A-U]“*F.  (27) 


When  the  coefficients  are  determined,  the  displacements  of  the  system  can  be  calculated  using 
equations  (7)  and  (8)  and  the  «  =  0  wave  propagation  coefficients. 
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The  elastic  plate  model  can  be  compared  and,  thus,  validated  for  a  thin  plate  at  low 
frequency  using  the  Bemoulli-Euler  differential  equation  of  motion  applied  to  a  thin,  infinite  flat 
plate  containing  periodic  masses.  This  model  has  been  previously  analyzed'^’  using  a  line  load 
excitation  and  is  reformulated  here  to  correspond  to  plane  wave  excitation.  This  expression  is 


=  -[Dk'^ -phco^)  '< 


1- 


Mco 


2  n=co 


z 


«  =  — 00 


H'l-l 


D 


k  + 


27m 
L  ) 


pho}"^ 


(28) 


where 


D  = 


Eh^ 

12(1 -f^^)  ’ 


(29) 


Figure  3  is  a  plot  of  the  transfer  function  of  displacement  in  the  z-direction  divided  by  input 
force  versus  wavenumber  at  a  frequency  of  200  Hz.  This  example  was  generated  with  the 
following  system  parameters:  thickness  is  0.01  meter,  density  is  1200  kg/m^.  Lame  constant  X  is 
2.25  X  10*  N/m^,  and  Lame  constant  p  is  2.50  x  lO’  N/m^.  In  figure  3,  the  dashed  line  is  thin  plate 
theory  without  the  periodic  masses  (i.e.,  M=  0);  the  solid  line  is  elastic  plate  theory  with  A/=  0.5 
kg/m  and  L  =  0. 1  meter  and  corresponds  to  equations  (8)  and  (27);  and  the  x  symbols  are  thin 
plate  theory  with  M=  0.5  kg/m  and  L  =  0. 1  meter  and  corresponds  to  equation  (28).  The  elastic 
plate  model  was  calculated  using  seven  modes  (-3  <  «  <  3)  that  produced  a  28-  x  28-element, 
system  matrix.  The  resonance  exhibited  by  the  thin  plate  model  is  the  plate  flexural  wave.  This 
energy  is  shifted  higher  in  wavenumber  as  masses  are  added  to  the  system.  Additionally,  the 
added  masses  allow  energy  in  the  system  at  characteristic  lengths — these  effects,  called  Floquet 
waves,  can  be  seen  at  around  20  rad/m  and  80  rad/m.  The  Floquet  waves  possess  stop-bands  and 
pass-bands  that  correspond  to  the  characteristic  lengths  of  the  structure.  Notice  in  figure  3  that  the 
thin  and  elastic  plate  models  with  periodic  masses  compare  closely,  but  not  identically.  This  is 
because  the  elastic  plate  equations  of  motion  allow  the  admittance  of  shear  wave  effects,  rotary 
inertia,  and  higher-order  wave  motion,  which  are  not  present  in  thin  plate  theory.  The  transfer 
function  in  the  jc-direction  is  not  shown  because  the  thin  plate  Bemoulli-Euler  equation  of  motion 
does  not  support  a  degree  of  freedom  in  this  longitudinal  direction. 
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Figure  3.  Transfer  Function  of  Displacement  Divided  by  Excitation  Force 

Versus  Wavenumber  at  200  Hz 


Figure  4  is  a  plot  of  the  elastic  plate  transfer  function  calculated  using  one,  three,  five,  and 
seven  terms  compared  to  thin  plate  theory  for  the  above  example.  The  wavenumber  axis  has 
been  expanded  to  250  rad/m.  This  plot  illustrates  the  behavior  of  the  model  as  more  terms  are 
added  to  the  truncated  model.  Notice  that  more  terms,  in  general,  correspond  to  dynamics  at 
higher  wavenumbers. 


Figure  4.  Comparison  of  the  Elastic  Plate  Transfer  Function  Calculated  Using  One,  Three, 
Five,  and  Seven  Terms  (Solid  Line)  and  the  Thin  Plate  Theory  (Dashed  Line) 
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3.  ELASTIC  PLATE  WITH  EMBEDDED  INTERIOR  MASSES 


The  elastic  plate  model  can  be  changed  so  that  the  periodic  masses  are  moved  from  the 
edge  into  the  interior  of  the  plate,  as  shown  in  figure  5. 


Figure  5.  Elastic  Plate  with  Periodic  Interior  Masses 


The  solution  to  this  problem  is  divided  into  two  parts  based  on  the  value  of  z.  The  upper 
region  of  the  plate,  location  b  <z  <c,is  denoted  with  a  prefix  superscript  (2),  and  the  lower 
region  of  the  plate,  location  a  <  z  <  b,  is  denoted  with  a  prefix  superscript  ( 1 ).  The  displacement 
in  the  x-direction  is 


u^{x,z,t) 


{k,  z,  co)exp(ikx)exp{io)t)  b  <  z  <c 
(k,z, co)exp{ikx)exp{icot)  a<  z  <b, 


where 


{k,  z,  co)  =  A{k,  co)ik  sin(Qrz)  +  B{k,  co)ik  cos(cirz) 
-  C{k,  co)P  cos(/?z)  +  D{k,  co)P  sin(y9z) 

and 


(30) 


(31) 


13 


^{k,z,o})  =  E{k,o})iks\n{az)  +  F{k,Q})ikcos,{az) 

-  G{k,co)/3cos{Pz)  +  H{k,co)Ps\r\{/3z). 

The  displacement  in  the  z-direction  is 

U,  {k,  z,  co)exp{ikx)exp{ia}t)  b<z<c 


u,{x,z,t)  =  > 


U,  (k,  z,  a>)exp(ikx)exp{ia)t)  a<  z  <  b. 


where 


^(k,z,a))  =  A(k,co)acos(az)  -  B(k,a))as\n(az) 

+  C(k,co)iksin{^z)  +  D(k,co)ikcos(/^z) 
and 


(k,  z,  co)  =  E{k,  co)a  cos(Qrz)  -  F{k,  co)a  sin(az) 

+  G{k,co)ik  sin(yffz)  +  H{k,co)ik  cos(/?z). 


(32) 


(33) 


(34) 


(35) 


There  are  eight  equations  of  motion  that  model  the  system.  The  first  two  equations,  i.e., 
(36)  and  (37),  are  the  normal  and  tangential  stress,  respectively,  at  the  top  surface  of  the  plate: 

^^h^^{k,c,co)  =  -F,  (36) 

and 

^'^^\^{k,c,co)  =  Q.  (37) 


The  next  four  are  the  interface  equations  at  the  plane  on  which  the  masses  reside:  two  (equations 
(38)  and  (39))  model  the  normal  and  tangential  stress  balances,  respectively,  between  the  upper 
and  lower  region  of  the  plate,  and  the  other  two  (equations  (40)  and  (41))  model  the  continuity  of 
displacement  in  the  x-  and  z-directions,  respectively.  These  interface  equations  are 
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=  Y"  ^{k +  ,b,co) , 

L  n=-oo  ^ 

(38) 

L  n=-oo  ^ 

(39) 

^'^'^U^ik,b,co)-^^^U^{k,b,co)^0, 

(40) 

and 

^{k,b,co)-^^'>U  ^{k,b,(o)  =  0. 

(41) 

The  final  two  equations  of  motion  are  the  normal  and  tangential  stress  at  the  bottom  surface  of 

the  plate: 

('^T^^(A:,a,ry)  =  0, 

(42) 

and 

(*>T^^(^,a,ry)  =  0. 

(43) 

Equations  (30)  through  (35)  are  now  substituted  into  equations  (36)  through  (43)  yielding 

[B(®\A:)]y(®^  =  [V^-*^(A:  +  2^)]y(±*)  +  g  , 

L 

(44) 

where  [B^°\^)]  is  an  8  X  8  matrix  that  models  the  dynamics  of  the  structure  for  «  =  0,  y^°^  is  the 

8  X  1  vector  of  wave  propagation  coefficients  for  «  =  0,  [V  (k  +  Inn  /  Z,)]  is  the  8  x  oo 

block- 

partitioned  matrix  that  represents  the  periodic  mass  loading  on  the  structure  for  «  =  -oo  to  «  =  qo, 
is  the  00  X  1  vector  of  wave  propagation  coefficients  for  n  =  -oo  to  «  =  oo,  and  g  is  the  8  x  1 
vector  that  models  the  plane  wave  excitation.  (The  matrix  and  vector  entries  in  equation  (44)  are 
listed  in  the  appendix.)  The  problem  solution  is  identical  to  that  shown  in  section  2  for  the  plate 
with  edge  masses,  although  each  A7-indexed  vector  has  eight  entries  associated  with  it  instead  of 

the  four  entries  when  the  masses  are  on  the  edge  of  the  plate.  Thus,  the  solution  becomes 

15 

y(±<»)  =[B-V]"‘G. 


(45) 


The  example  problem  of  an  elastic  plate  with  discrete  masses  at  the  bottom  edge  was  re¬ 
analyzed  and  the  displacement  results  were  identical  in  the  z-direction  and  nearly  identical  in  the 
x-direction.  This  indicates  that  moving  the  masses  into  the  interior  has  no  effect  on  the  z- 
displacement  and  only  a  slight  effect  on  the  x-displacement  at  this  frequency,  which  is  an 
expected  result  because  this  is  a  low-frequency  example  where  the  wavelengths  of  the  forcing 
function  are  much  larger  than  the  thickness  of  the  plate. 
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4.  NUMERICAL  EXAMPLE 


A  numerical  example  is  now  analyzed  to  show  the  dynamic  response  of  elastic  plates  with 
masses  on  the  edge  and  in  the  interior.  This  example  was  generated  with  the  following  system 
parameters:  thickness  is  0. 1  meter,  density  is  1200  kg/m\  Lame  constant  A  is  4.5  x  10®  nW, 
Lame  constant  /y  is  5.0  x  10^  N/m^,  mass  per  unit  length  is  0.5  kg/m,  and  the  mass  separation 
distance  is  0.3  meter.  Both  elastic  plate  models  were  calculated  using  15  modes  (-7  £  «  <  7)  that 
produced  a  60-  x  60-element  system  matrix  for  the  edge  mass  problem  and  a  120-  x  120-element 
system  matrix  for  the  interior  mass  problem.  For  the  interior  mass  problem,  the  masses  were 
located  at  the  mid-plane  of  the  plate.  The  displacement  values  were  output  at  location  of  0.25h 
(0.025  meter)  from  the  top  of  the  plate.  This  problem  is  constructed  so  that  the  wavelengths  of 
the  shear  and  dilatational  waves  (0.204  and  0.677  meter,  respectively)  were  on  the  order  of  the 
length  scales  of  the  plate  and  mass  separation  distance. 

Figure  6  is  a  plot  of  the  dispersion  curve  for  the  plate  with  edge  masses,  which  corresponds 
to  free-wave  propagation  of  the  system.  A  frequency  range  of  0-2000  Hz  was  chosen  so  that  the 
first  three  waves  are  present  in  the  plot.  The  three  dark,  solid  lines  represent  locations  in  the 
wavenumber-frequency  plane  waves  that  are  centered  about  ^  =  0.  These  waves  are  also  present 
in  the  infinite  plate  without  periodic  masses.  The  other  waves  (Floquet  waves),  shown  as  dashed 
lines,  correspond  to  waves  that  share  an  integer  multiple  of  a  characteristic  wavenumber.  This 
wavenumber  is  determined  by 

,  In' 

kc=—  (46) 

and  is  equal  to  20.9  rad/m.  These  Floquet  waves  represent  energy  that  is  (partially)  trapped  in 
the  spaces  between  the  periodic  masses  and  is  converted  into  additional  free-wave  propagation. 

Figures  7  and  8  present  the  transfer  function  of  tangential  displacement  and  normal 
displacement,  respectively,  divided  by  excitation  force  versus  frequency  and  wavenumber  for  the 
plate  with  edge  masses.  In  both  figures,  the  data  are  displayed  in  the  decibel  scale  referenced  to 
meters  per  Pascal.  The  free  wave  that  appears  in  the  dispersion  curve  is  clearly  evident  in  both 
the  tangential  (figure  7)  and  normal  (figure  8)  displacement  data.  Figures  9  and  10  present  the 
transfer  function  between  tangential  displacement  and  the  normal  displacement,  respectively,  and 
force  versus  wavenumber  at  1000  Hz.  Figures  9  and  10  are  divided  into  two  plots  for  clarity;  the 
upper  plot  is  the  problem  of  an  elastic  plate  with  edge  masses  and  the  lower  plot  is  the  elastic 
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plate  with  interior  masses.  The  dashed  line  in  all  four  plots  is  the  elastic  plate  solution  with  the 
absence  of  masses  (A/=  0). 

Upon  examination  of  figures  9  and  10,  several  features  are  noted.  The  energy  field  has 
become  modally  dense  as  the  mass  (with  its  characteristic  length)  now  introduces  additional 
resonances  and  anti-resonances  in  the  wavenumber  space.  This  result  is  expected  because  the 
model  has  changed  from  a  continuous,  homogeneous  infinite  structure  to  a  structure  that 
periodically  inhibits  wave  motion  by  the  addition  of  discrete  masses.  Additionally,  this  effect 
has  also  been  evident  in  previously  developed  low-frequency  (bending-wave)  models.  The 
overall  energy  levels  are  approximately  the  same  between  the  models  with  masses  and  the 
models  without  masses.  This  is  because  higher-frequency  models  are  influenced  more  by 
stiffness  of  the  system  than  by  mass.  Mass  effects  are  seen  more  predominately  at  lower 
frequencies  (see  figure  3). 
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Figure  6.  Dispersion  Curve  of  Plate  with  Edge  Masses 
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Figure  7.  Transfer  Function  of  Tangential  Displacement  Divided  by  Excitation  Force  Versus 
Frequency  and  Wavenumber  for  a  Plate  with  Edge  Masses 
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Figure  8.  Transfer  Function  of  Normal  Displacement  Divided  by  Excitation  Force  Versus 
Frequency  and  Wavenumber  for  a  Plate  with  Edge  Masses 
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Figure  9.  Transfer  Function  of  Tangential  Displacement  Divided  by  Excitation  Force  Versus 
Wavenumber  at  1000  Hz  for  Edge  Masses  (Upper  Plot)  and  Interior  Masses  (Lower  Plot) 
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Figure  10.  Transfer  Function  of  Normal  Displacement  Divided  by  Excitation  Force  Versus 
Wavenumber  at  1000  Hz  for  Edge  Masses  (Upper  Plot)  and  Interior  Masses  (Lower  Plot) 
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5.  CONCLUSIONS 


The  solution  of  an  elastic  plate  containing  periodic  edge  and  embedded  masses  has  been 
derived  and  found  to  compare  favorably  with  previously  developed  thin  plate  models  at  low 
frequency.  A  numerical  example  of  high-frequency  dynamics  was  presented  and  the  details 
discussed.  The  dispersion  curve  and  transfer  functions  of  tangential  and  normal  displacements 
were  illustrated.  It  was  shown  that  the  lower-  and  higher-frequency  waves  propagate  at  spatial 
lengths  that  correspond  to  integer  multiples  of  the  separation  distance  of  the  periodic  masses. 
This  characteristic  makes  the  system  become  modally  dense,  even  at  low  frequency  and  low 
wavenumber. 
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APPENDIX 

MATRIX  AND  VECTOR  ENTRIES 


The  entries  of  the  matrixes  in  equations  (20)  and  (21)  are  listed  below.  Without  loss  of 
generality,  the  bottom  of  the  plate  is  defined  as  z  =  o  =  0.  For  the  +  27m ! L)]  matrix,  the 

nonzero  entries  are 

"11  =(-«^^^-2«^-"-^^;J)sin(«^/;)  , 

(A-1) 

"12  =(-««^-2«^y"-^^)cos(a„/7)  , 

(A-2) 

"13  =  cosifi^h)  , 

(A-3) 

«14  =  sin(y5„/7)  , 

(A-4) 

oji  =2i/jk„a^cos{a^h)  , 

(A-5) 

^22  =  S>n(«„/?)  , 

(A-6) 

"23  -^^)sin(/9„/i)  , 

(A-7) 

^24  =  cosi/3 ^h)  , 

(A-8) 

"32  . 

(A-9) 

"33  =  , 

(A- 10) 

A-1 


«4i  . 

(A-11) 

and 

^44  ’ 

(A- 12) 

where 

(A- 13) 

(A- 14) 

and 

(A- 15) 

For  the  [U^”^(A:  +  2nnl  L)]  matrix,  the  nonzero  entries  are 

-Mo? 

"31  =  1 

(A- 16) 

-Mo? 

"34  =  ^  ’ 

(A- 17) 

"42  =  ^  ’ 

(A- 18) 

and 

(A- 19) 

A-2 

The 


vector  entries  are 


x^®^  =[/l(A:,<y)  B{k,(o)  C{k,co)  D{k,(o)Y 


^(0)  ^(0)  ^(0)  ^(0) 


(A-20) 


The  x^-*) 


vector  entries  are 


T±«>)  _ 


=  [... 


(-1)  ^(-l)  r(-l)  n(-l)  ^(0)  ft(0) 


C 


c(0)  £>(•>  ...j''. 


(A-21) 


The  f  vector  entries  are 

f  =  [-F  0  0  0]^.  (A-22) 

The  entries  of  the  matrixes  in  equation  (44)  are  listed  below.  Without  loss  of  generality, 
the  location  of  the  masses  in  the  z-direction  is  defined  as  z  =  6  =  0.  For  the  +  Ittki!  L)] 

matrix,  the  nonzero  entries  are 


^11  =(-«^^-2Qr^y“-/l^^)sin(a„c)  , 

(A-23) 

^12  =(-«^^-2«^/^-/W:^)cos(a„c)  , 

(A-24) 

^13  =  cos(/?^c)  , 

(A-25) 

^14  =  sin(;0^c)  , 

(A-26) 

^21  =  cos(or^c)  , 

(A-27) 

^22  =  sin(a^c)  , 

(A-28) 

A-3 


^23  -kl)sxx\{p^c)  , 

(A-29) 

^24  -kl)cos{P^c)  , 

(A-30) 

(A-31) 

hi  =  , 

(A-32) 

he  +  +  ’ 

(A-33) 

^37  =  ’ 

(A-34) 

Z?4i  =  2\/jk^a^  , 

(A-35) 

'>44=MA„^-*^, 

(A-36) 

^5  =  , 

(A-37) 

*48=-M/?„'-^')  . 

(A-38) 

II 

-C 

(A-39) 

^54  =  ’ 

(A-40) 

s: 

1 

II 

-C 

(A-41) 

^58  =  > 

(A-42) 

A-4 


^62  =  > 

(A-43) 

^63  =  ~f^n  ’ 

(A-44) 

^66  =  ’ 

(A-45) 

^67  =  ’ 

(A-46) 

^>75  =(-Qr^A-2Qr^/^-/W:^)sin(or^<3)  , 

(A-47) 

^76  =(-««^-2or^y“-^^^)cos(a„o)  , 

(A-48) 

b-n  =  ^i^nPn  cos(/0„a)  , 

(A-49) 

^78  =  sin(y9^o)  , 

(A-50) 

^85  =  COS(«„a)  , 

(A-51) 

^86  =  sin(«„a)  , 

(A-52) 

^87  -^^)sin(y0„a)  , 

(A-53) 

and 

b%%=M{Pl  -kl)cos{p^a)  . 

(A-54) 

The  +  2nnl  L)]  matrix  can  be  written  as 

A-5 

(A-55) 


2m 


)  0 


where  +  Iren! L)]  is  defined  in  equations  (A-16)  through  (A-19)  and  0  is  a  4  x  4  matrix 

whose  entries  are  all  zero.  The  vector  entries  are 


=[A{k,(o)  B{k,co)  C{k,(o)  D{k,co)  E{k,co)  F{k,co)  G{k,co)  H{k,(o)Y 


^(0)  ^(0)  ^(0)  ^(0)  ^(0)  ^(0)  (j(0)  ^(0) 


(A-56) 


The  vector  entries  are 


,(±«>)  _ 


^(-1)  ^(-1)  (^(-1)  £)(-l)  £(-l)  p{-\)  ^(-1) 

^(0)  ^(0)  ^(0)  ^(0)  ^(0)  ^(0)  ^(0)  ^(0) 

^(1)  5(1)  d(*>  //(•)  - 


(A-57) 


The  g  vector  entries  are 


g=[-F  0  0  0  0  0  0  0]  . 


(A-58) 


A-6 
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